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O I Abstract. We study the factorised steady state of a general class of mass transport 

models in which mass, a conserved quantity, is transferred stochastically between sites. 
Condensation in such models is exhibited when above a critical mass density the marginal 
distribution for the mass at a single site develops a bump, p C ond(™), at large mass m. This 
bump corresponds to a condensate site carrying a finite fraction of the mass in the system. 
Here, we study the condensation transition from a different aspect, that of extreme value 
statistics. We consider the cumulative distribution of the largest mass in the system and 
compute its asymptotic behaviour. We show 3 distinct behaviours: at subcritical densities 
the distribution is Gumbel; at the critical density the distribution is Frcchet, and above 

ch ■ the critical density a different distribution emerges. We relate p CO nd(rn) to the probability 

O . density of the largest mass in the system. 
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PACS numbers: 05.70.Fh, 02.50.Ey, 64.60.-i 



1. Introduction 

Mass transport models form a general class of nonequilibrium systems where some conserved 
quantity, which we will refer to as mass, is transported stochastically between the sites of a 
lattice, according to certain prescribed dynamical rules [1] . At various levels of description 
such models may represent the dynamics of traffic flow [2], [3], granular clustering [4], phase 
ordering [5], network rewiring [H1EIE], force propagation [9], aggregation and fragmentation 
[TUt [TT] and energy transport [12] — for a review see [13] . 

Of particular interest is the nonequilibrium steady state which is attained in the long 
time limit. In general the structure of nonequilibrium steady states is not known, however 
in some cases the steady state factorises. That is, the joint distribution of the collection of 
masses {mj at sites i = 1, . . . ,L is given by 

where Z(M, L) is just the normalization 

L poo I L \ 

Z(M,L) = H dmJK^mpM (2) 

and is the analogue of the 'canonical partition function'. Note that (pQ) is a product of 
single-site weights f{rrii) but the 5-function in ([I]) imposes the global constraint that the 
total mass, M, in the system is conserved: 

M = 5>,. (3) 

i=i 

Thus, correlations are induced between sites and in general the single-site mass probability 
distribution, i.e., the marginal p(m, M, L) = J P(m, mi-, m^, . . . , mi) Y\f =2 drrii =£ f(m). (For 
brevity we will often use p(m) as a shorthand for p(m, M, L).) 

In the case of a simple class of one- dimensional asymmetric mass transport models, a 
necessary and sufficient condition for factorisation has been determined. In these models the 
dynamics is defined as follows. At each time step, a portion, rhi < m% of the mass at each 
site, is chosen from a distribution (f)(jhi\mi) and is transferred to site % + 1. The stationary 
state is factorised provided that the kernel (p(rh\m) is of the form 

, . . . u(fh)v(m — m) 
0(m|m) = f(m) ' ( } 

where u(z) and v(z) are arbitrary non-negative functions. The single-site weight is then 
given by 

f(m) = / dm u(rh)v(m — m) . (5) 

Jo 



This model is general enough to include many well-known models as special cases pQ. 
Choosing the chipping kernel 0(m|m) appropriately, recovers the Zero-Range Process 
(discrete, integer valued masses with m restricted to 0,1) which always has a factorised 
steady state, and the Asymmetric Random Average Process (continuous masses with m 
chosen as a random fraction of m) which exhibits a factorised steady state in certain cases 
[HI [HI [15] . Moreover, the model encompasses both discrete and continuous time dynamics. 
The condition for factorisation has been generalised to arbitrary dimensions and arbitrary 
graphs [16] . where generally one is able to demonstrate sufficient conditions for factorisation. 
We also note that factorised steady states have been extensively studied in the queueing 
theory literature, see e.g. [17] . 

By choosing the dynamics appropriately various forms for the steady state weight f(m) 
can be produced from ([5]). Of particular interest has been the case where for large m 

f(m) ~ AmT 1 . (6) 

In the following we choose the normalisation constant A so that f£° dm f(m) = 1. Then, if 
the index 7 > 2 the phenomenon of condensation occurs ^8] l2l [T%] [20l [2T| [22l [23| IM| [25]. 
This is manifested by three distinct large m behaviours of p(m) as the global mass density, 
p, given by 

M m 

P= T , (7) 

is increased (see Figured]): 

(i) When p < p c , p(m) ~ Bf(m)z m where z < 1 and B is a normalisation constant for 
p(m). This is the fluid phase where p(m) decays exponentially for large enough m 

(ii) When p = p c , p(m) ~ f{in) — Am' 1 . This is the critical point where p(m) decays as a 
power law. The critical density p c is given by 

p c = dmmfOm) . (8) 

JO 

(iii) When p > p c , for finite but large L, p{m) develops an additional bump centred at 
m = L(p — p c ). Roughly speaking, for large L, p{m) ~ f(m) + p C ond(^, M, L). This is 
the condensed phase where this extra piece of p(m), p C ond( m , M, L), emerges in addition 
to the critical point distribution. This piece (the bump in p(m)) carries a weight 1/L. 
It thus represents a single site containing the excess mass L(p — p c ), which coexists with 
a background of power-law distributed critical mass. The large system-size behaviour 
of the bump p C ond(^, M, L) was determined [23], El] and shown to have distinct forms 
in the regimes 7 > 3, where on the scale m — L(p — p c ) = 0(L 1//2 ) it is gaussian, and 
2 < 7 < 3, where the distribution is highly asymmetric and non-gaussian. These two 
forms imply that the number of particles in the condensate has fluctuations of 0(L 1 / 2 ) 
for 7 > 3 but anomalously large fluctuations of (^(L 1 ^ 7-1 )) for 3 > 7 > 2. 



-2 



Q. 



-7 - 




-12 



L=100 



1 3 

ln(m) 



p= 

I 

5 



Figure 1. The exact single-site mass distribution p(m) plotted for a particular choice of 
f(m) with 7 = 5/2 and p c — 2, and system size L — 100: full line p = 1 (subcritical); dotted 
line p = 2 (critical); dashed line p = 6 (supercritical). The condensate bump, Pcond( m )i is 
apparent in the supercitical curve. 



This description of condensation has been from the point of view of the marginal distribution 
p(m) [231 [2l] i n the thermodynamic (large L) limit. However, condensation is a co-operative 
behaviour resulting in the emergence of a single extensive mass amongst the L masses. While 
the signature of the condensation transition is manifest in the single-site mass distribution 
p{m), a more direct and somewhat natural probe for the condensation is the study of the 
statistics of the maximal mass in the system. This is because in the condensed phase, the 
site where condensation occurs is also the site carrying the maximal mass. It is from this 
perspective, that of the extreme value statistics of the masses, that we study condensation 
in this work. 



2. Condensation and extreme value statistics 

In this paper we ask the question — what is the distribution of the largest mass in the system? 
Thus rather than the single-site mass distribution p(jn, M, L) studied previously, we consider 
the probability Q(jn, M, L) that the largest mass in the system (of size L and total mass M) 
is less than or equal to m, which is an example of an extreme value distribution. 

In the zero-range process the asymptotic size of the largest mass m* has been considered 
by Jeon et al [20]. In that work the case ([6]) was considered with 7 > 3 and it was shown 
that in the fluid phase m* = O(lnL), at the critical density m* = 0(L 1 ^' y ~ 1 ' > ), and in the 
condensed phase m* / L = p — p c . Here we study the statistics of the maximal mass in the 



more general mass transport model described in the introduction. Moreover, we obtain the 
full asymptotic probability distribution of the maximal mass for large L and not just its 
average size. 

The theory of extreme value statistics (EVS) is well understood for the case of 
independent random variables. If one considers a set of L independent and identically 
distributed (i.i.d) random variables, each drawn from a common distribution f(m) then the 
limiting distribution of their maximum, appropriately centred and scaled, has three possible 
forms [26]. The limiting distribution is: a Gumbel distribution when m unbounded and f(m) 
decays faster than a power law for large m; a Frechet distribution when m is unbounded and 
f(m) decays as a power law for large m, and a Weibull distribution when m is bounded. 

In our case, however, the L masses rrii are correlated due to the global constraint of 
conserved total mass, explicitly manifest in the delta function in (Tj[|). Had this delta function 
constraint been not there, since f(m) ~ Am~ y , one would expect a Frechet distribution for 
the maximal mass based on the EVS of i.i.d. random variables. Here, we find that the delta 
function induces nontrivial correlations that modify this naive expectation based on the EVS 
of i.i.d. random variables. 

We shall show that the three phases (respectively for p < p c , p = p c and p > p c ) 
exhibit distinct forms for Q(m, M, L), the cumulative probability distribution for the largest 
mass. In the fluid phase (p < p c ) Q(m, M, L), appropriately centred and scaled, is a Gumbel 
distribution; at the critical point (p = p c ) it becomes a Frechet distribution; in the condensed 
phase (p > p c ) the distribution of the largest mass in the system is given by Lp conc i(m, M, L) 
where p con ^(m,M,L) was computed explicitly in [23], [21]. m the latter case we see that 
maximal mass distribution is essentially given by the tail of the single-site mass distribution. 
The multiplying factor L simply denotes the fact that any of the L sites can carry the 
maximal mass. Thus, in the condensed phase, the maximal mass distribution is neither 
of the three types (Gumbel, Frechet or Weibull) that one encounters in the EVS of i.i.d 
random variables, indicating strong effects of the correlations between the masses. Thus the 
correlation arising from the global mass conservation modifies the maximal mass distribution 
both in the fluid phase (p < p c ) as well as in the condensed phase (p > p c ). Interestingly 
exactly at the critical point p = p c , the maximal mass distribution is Frechet as would be 
prediction based on i.i.d variables indicating that the correlation is somehow least effective 
exactly at the critical point. 

3. Computation of the distribution of the largest mass 

We begin by defining the cumulative distribution Q(x, M, L), the probability that the largest 
mass is less than or equal to x: 



Here, 



I(x, M,L) = f[ f dmj(m t ) 5(£ m, - M) 

i=\ J ° 7=1 



(10) 



and Z(M, L) (gj) is given by Z(M, L) = J(oo, M, L). Equation © follows from the fact that 
if the largest mass is < x, all masses must be < x. We shall also consider the probability 
density of the largest mass 



P(x,M,L) = ^-Q(x,M,L) 



Lf(x 



I(x,M -x,L- 1) 



Z(M,L) 

Expression (Tl2l may be compared to the single site mass distribution studied in 

, w rN ,, ,Z(M-m,L-l) 

p(m,M,L) = f(m)— ; — —, . 

^ v ' ' ; M J Z{M,L) 

The Laplace transform of I(x, M, L) is easy to compute from fllOp 



dM I(x,M,L)e 



-sM 



dm f(m)e 



therefore 



I{x, M, L) 



c+oo £s 



exp L I ps + In 



dm f(m)e 



(11) 
(12) 

1J 
(13) 

(14) 
(15) 



/c-joo 2rci 

where M = pL and c is chosen so that integration contour is to the right of any singularity 
of the integrand. Our task now is to evaluate the integral ffT5|) . In [231 1211 the partition 
function Z(M, L) = I(oo, M, L) was evaluated and we shall follow a similar approach here. 



3.1. Fluid phase p < p c 

First we consider the fluid phase in which case I(x, M, L) may be evaluated by the saddle- 
point method. Let us define 



H"n\Xy S ) 



dm f(m) m n e 



(16) 



where due to our choice of normalisation in (jBJ), yu (0, 0) = 1. 
For < s < oo and x large, p n (x, s) is asymptotically 

p n ( x ,s) = J -^ [1 + 0(1/*)]. 

s 

On the other hand, for s = and f(m) ~ Am -1 for large m, we have for large x, 

A 



(17) 



p n (x, 0) 



-x 



n ~ 7+1 for 7>n + l. 



7 — n — 1 

Next, using f£ dmf(m)e~ sm = /Xo(0, s) — /io(^, s ), we ma Y write (TT5"j) . for large x, as 

/i (x, s 



:u 



J(x, M, L) ~ / — — expL < ps + ln/x (0, s) 

Jc-ioD 2m [ 



A*o(0,s' 



+ 



(19) 
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and to leading order the saddle point, s , of the integrand is independent of x and is given 
by the equation 

jUi(0,s ) 



P 



/i (0, s ) 
d 



(20) 



where we have used — [/io(0, s)l = — /ii(0, s). The saddle point exists and lies on the positive 

as 
real axis when p < p c given from ([8]) by 

p c = Pi(0,s). (21) 

Therefore, in the fluid phase (p < p c ), to leading order, we obtain 



Q(x,M,L) ~ exp 



Po{x,s ) 
Vo(0, so) 



exp < — L 



Ae~ 



x^sgiso) J 



(22) 



where g(so) = /io(0, so), and we have used (TT71) . For x = — — hO(lnlnL) this may be written 
in the familiar form of a Gumbel distribution 

' x — a>L 



Q(x, M, L) ~ exp < — exp 



where 



ol = — In 

so 



LA 



1/so 



(23) 



(24) 



.(lnL)7 S ^( So ) 

Thus the probability density of the maximal mass P(x, M, L) has a peak around x = a^ 
with a width 6^. The scaling form of the distribution around this peak has the Gumbel form 
and implies that the largest mass is O(lnL). 



3. 2. Critical density p = p c 

At the critical density p = p c the saddle point of the integral is at so = and p c is given by 
(JH]). In this case, we use (TTS1) 

A 1 



p (x, 0) 



7- laf" 1 



.4 



and find, to leading order, 

Q(x,M,L)^exp|-^ (o)(7 _ i) j 

This is a Frechet distribution and implies that the largest mass is 0(L 1 / < - 7-1 ^) 



-Lx-^'-V 



(25) 



(26) 



3.3. Condensed phase p > p c 

We now turn to the condensed phase where p > p c . In this case the saddle-point method 
can no longer be used to evaluate f|T5|) as there is no longer a solution to fl20l) . 

It turns out to be convenient, in this case, to work directly with the probability density 
of the largest mass P(x,M,L) in (Tl2|) . rather than the cumulative distribution Q(x,M,L). 
In fact, the calculation of ( fT2l) reduces to precisely that of the distribution p con d(m,M,L) 
detailed in [231 121] which we now review. 

To analyse ( TT2l) . we need to consider both I(x,M — x, L— 1) (in the numerator) and 
also Z(M,L) (in the denominator). 

First, the behaviour of Z(M, L) = /(oo, M, L) for p > p c has been computed in [23|l2l]. 
The integral (TTo^) 

— -expL[ps + ln/i (0,s)] (27) 

-ioo 2m 

will be dominated by small values of s. Therefore, one expands the term ln/z (0, s) for s 
small. The expansion depends on the value of 7 and is given by 

For 3<7, ln/x (0,s) ~ - sp c + bs^ 1 + ■ ■ ■ (28) 

For 2<7<3, ln/x (0, s) ~ - sp c + — s 2 + ■ ■ ■ + bs"'' 1 + ■ ■ ■ (29) 



w „, t- ^ f io ° ds exi)LH(s,x) , . 

I(x,M -x,L-l)= / , 30 

v ' ' ; J-ioo2m u n 0,s v ; 



where b = AT(1 - 7), A 2 = (/x 2 (0, 0) - /i 2 (0, 0))//i (0, 0) and p c is given by (G5J1. 

Similarly, we consider the numerator in ffT2l) . for which the relevant integral ffl5l) is 

r io ° ds expLiJ(s,x) 
/x (0,s) 
where for large x (we will be considering x = O(L) in the large L limit) 

H(s,x) ~ ps — xs/L + ln/i (0, s) . (31) 

We now expand this expression for H(s,x) for small s, using (I28II29I) . 

As a result of the small large L, small s expansions just described, one obtains from 

P^M.L^L/M ^'ff-^-f (32) 

W L {-{p- p c )) 

where Wi,{y) has different forms according to the cases 7 > 3 and 2 < 7 < 3 which we 
discuss separately below. 

We note that fl32l) yields the same leading large x behaviour for P(x,M,L)/L as for 
the single site mass distribution p(x, M, L) computed in [21] . This is due to the fact 
that for 1 > lni the same leading behaviour is obtained for I(x, M — x, L — 1) as for 
I(oo, M — x, L — 1). Moreover, we can identify the peak in the probability density of the 
largest mass with the peak p CO nd in the single site distribution. 



Case (i) 7 > 3 

In this case Wl(u) in (1521) is given by 

/too Ac r A 

— expL - ys + - s 2 /2 + --- + fo^ 1 + -- 
-ioo ^7T2 I 



(33) 



The important behaviour of the function Wl{d) determined in [23l [24] may be summarised 
as 



W L (y) 
W L (y) 



1 



v^ttA 2 !, 
A 



-y 2 L/2A 



ly^Li- 1 
from which we deduce that for p > p c 



for y = o(L- 1/3 ) 
for y = 0(l), 



W L {x/L-{p-p c )) 



(34) 
(35) 

(36) 



Thus the probability density of the maximal mass is peaked around x = (p — p c )L and near 
this peak, that is, over a scale of x — L(p — p c ) = o(L 2//3 ), the probability density has a 
normalized gaussian form 



P(x,M = pL,L) 



V2nA 2 L 



cxp 



(x-(p- p c )L)' 
2A 2 L 



(37) 



Case (ii) 3 > 7 > 2 

From (fT2l one obtains (I3"2l where now 

W L (y) = L- 1 ^-^V y [L^- 2 ^- 1 ^y 

and 



Yr(z) 



ioo ds 



(38) 



(39) 



-»oo 27TZ 

with, as before, b = AT(1 — 7). The function V^(z) is highly asymmetric around z = and 
has the following asymptotic behaviour 

V 7 (*) * ^ 



as 2; — > — 00 



as 2 — > 00 



~ ClZ (3~7)/2(7-2) expf-caz^ 1 )/^- 2 )] 

where ci and c 2 are two constants that were computed in 

d= [27r(7-2)(6( 7 -l)) 1/( ^ 2) ] _1/2 , 

c 2 = (7-2)[(7-l)(6(7-l)) 1/(7 - 2) ] _1 • 
Using (1381 and (T4"0"j) in (I3"2"j) . we obtain the normalized maximal mass density 



P(x, M, L) 



£1/(7-1) V t 



x - L(p - p c 

L 1 /(7-i) 

9 



(40) 

(41) 

(42) 
(43) 

(44) 



Thus, in the condensed phase (p > p c ) the maximal mass distribution has neither of the three 
forms associated with the EVS of i.i.d random variables. Rather, the scaled distribution has 
a gaussian form for 7 > 3 and has a nontrivial non-gaussian form for 2 < 7 < 3. In both 
cases the largest mass is O(L). 

4. Discussion 

In summary, we have probed the condensation transition in a generalized class of mass 
transport models by studying the cumulative probability distribution of the maximal mass in 
the system. The masses in this model are globally constrained via the mass conservation law. 
Without this global constraint the masses would have been independent random variables 
each drawn from a power-law distribution and based on the EVS of the i.i.d random variables 
one would expect a Frechet distribution for the maximal mass for all densities p. Instead, 
we have shown via exact asymptotic calculation that the global constraint is sufficiently 
strong to modify this expectation based on i.i.d. variables and one obtains respectively a 
Gumbel distribution (for p < p c ), a Frechet distribution (for p = p c ) and a completely 
different distribution in the condensed phase (for p > p c ). In the latter case, the distribution 
is gaussian for 7 > 3 and highly non-gaussian for 2 < 7 < 3. Thus, our exact results for 
this model form a useful addition to the list of exactly solvable cases for the distribution 
of extreme values of a set of correlated random variables, for example, the eigenvalues of 
random matrices [271 1211 [291 USB EH [32] , energies of configurations a directed polymer in a 
random medium [33] , the heights in a one-dimensional interface [3U [35], [361 EZ] • We note 
that the Gumbel distribution emerging out of a global constraint was also observed recently 
in the context of complex chaotic states [3H] • 

We have also seen how in the condensed phase the bump, p con d( m ) in the single-site 
mass distribution is directly related to the probability density of the size of largest mass in 
the system. 

Finally, we note that as the density p is increased the phenomenon of condensation is 
manifested through changes in the asymptotic behaviour of the large mass distribution as 
described above. It would be of interest to analyse more closely the crossovers, both from 
the Gumbel distribution for p < p c to the Frechet distribution at p = p c and from the 
Frechet distribution at p = p c to the condensed phase distribution for p > p c , as the density 
is increased. This would involve the consideration of finite-size effects and finite-size scaling 
near the transition. 
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